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\^Summary 

The  rapid  crack  propagation,  crack  curving  and  arrest  mechanisms  associated  with 
a  pressurized,  thin*walied  ductile  steel  tubes  are  used  to  develop  a  model  of  axial 
rupture  of  an  aircraft  fuselage.  This  model  is  used  to  replicate  axial  crack 
propagation  along  a  line  of  multi-site  damage  (MSD)  and  crack  curving  and  arrest 
near  a  tear  strap  of  an  idealized  fuselage.  ^ 

Introduction 

"Static  linear  elastic  fracture  mechanics  (LEFM)  has  been  used  to  predict  a 
monotonically  increasing  stress  intensity  factor  with  increasing  crack  iength  of 
an  axial  crack  in  a  pressurized  pipe  [1,2]  as  well  as  to  model  dynamic  crack 
branching  in  a  pressurized  steel  pipe  [3].  Plasticity  effect  has  also  been 
incorporated  to  this  LEFM  analysis  using  semi-empirical  adjustments  [4,5],  or  by 
kinematic  modelling  of  the  yielded  pipe  [6].  Many  of  these  analyses,  however,  do 
not  consider  the  dynamic  effects  generated  by  the  propagating  axial  crack  as  well 
as  the  influence  of  the  large  plastic  deformation  of  the  crack  flaps. 

Under  the  sponsorship  of  the  US  Department  of  Transportation  (DOT),  one  of  the 
author  and  his  colleagues  studied  axial  crack  propagation,  curving  and  arrest  in  a 
small  scale  model  of  a  line  pipe  [7,8].  These  experimental  results  were  then  used 
to  modify  a  numerical  model,  which  was  previously  developed  for  analyzing 
bursting  thick-wailed  steam  pipes.  This  study  showed  that  the  main  driving  force 
of  the  axial  crack  in  a  bursting  pipe  was  the  gas  pressure  acting  on  the  crack  flaps 
which  literally  tore  the  pipe  apart.  It  was  also  found  that  the  axial  stretch  of  the 
crack  flap,  not  only  served  as  an  energy  dissipation  mechanism,  but  also  imposed 
an  axial  tensile  stress,  which  is  larger  than  the  hoop  stress,  as  much  as  five  pipe 
diameters  ahead  of  the  crack  tip  thus  leading  to  crack  curving  and  crack  arrest. 

While  an  aircraft  fuselage  is  an  infiniteiy  more  complicated  structure  than  a 
monolithic  line  pipe,  the  crack  propagation,  curving  and  arrest  mechanisms  of  the 
latter  could  be  of  use  in  a  postmortem  analysis  of  a  fuselage  failure  as  well  as  for 
estimating  the  crack  arrest  capability  of  a  tear  strap,  in  the  following,  a  brief 


1 


review  of  the  cogent  results  generated  in  the  DOT  pipe  fracture  study  as  well  as 
the  crack  curving  and  crack  branching  criteria,  which  were  developed  under  a 
Office  Naval  Research  (ONR),  contract  are  presented.  These  results  are  then  used 
to  assess  the  crack  arrest  capability  of  a  hypothetical  fuselage  with  tear  straps. 

Analysis  of  Axial  Rupture  of  Subsize  Gas  Transmission  Lines 

The  experimental  phase  of  the  DOT  studies  consisted  of  measuring  the  dynamic 
depressurization  rate,  crack  extension  rate,  static  and  the  dynamic  axial  and 
circumferential  strains  along  the  axial  crack  in  a  2-in.  diameter,  schedule  10, 
carbon  steel  pipe.  These  pipes  were  pressurized  to  approximately  80  percent  of 
the  yield  strength  in  hoop  stress  in  a  series  of  static  and  dynamic  pipe  rupture 
experiments.  The  test  section  ‘of  the  pipe  was  pre-grooved  to  a  depth  of 
approximately  20  percent  of  the  pipe  wall  thickness  to  simulate  a  locally 
embrittled  zone  for  continuous  and  controlled  axial  crack  propagation  in  the 
ductile  pipe  wall.  Such  grooving  simulated  multi-site  damage  (MSO)  in  a  fuselage. 
The  grooves  were  carefully  machined  to  avoid  local  elevation  of  the  ultimate 
strength  of  the  pipe  wail  material  due  to  strain  hardening,  and  thus  the  pipe  wall 
fractured  with  a  100%  shear  lip.  The  propagating  crack  always  arrested  upon 
entering  the  ungrooved  section,  which  simulated  the  fuselage  section  with  a  tear 
strap  of  the  pressurized  pipe.  The  rupturing  profiles,  i.e.,  the  pipe  flap  motions, 
were  recorded  by  a  higfr  speed  framing  camera.  Rgure  1  shows  typical  sequential 
photographs  of  a  rupturing  pipe  where  the  crack  bifucated  upon  entering  the 
ungrooved  section  of  the  pipe. 


Fig.  1.  Axial  crack  bifucation  and  arrest  in  a  pressurized  steel  pipe 

A  special  software  for  processing  the  high-speed  photographic  records  and  for 
determining  the  crack  opening  displacements  (COD)  and  axial  straining  along  the 
crack  flaps  was  developed.  This  data  was  used  to  compute  the  crack  tip  opening 
angle  (CTOA).  The  static  and  dynamic  crack  opening  shapes  and  CTOA's  were 
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similar  and  the  CTOA  remained  constant  over  a  crack  velocity  range  of  600  to  1000 
fps  but  varied  with  the  groove  depth.  When  the  groove  depth  was  decreased  to 
simulate  a  tougher  pipe  material,  the  CTOA  increased  and  the  crack  velocity 
decreased. 

Records  of  the  strain  gage  rosettes,  which  were  placed  along  the  crack  path  in  the 
dynamic  pipe  rupture  tests,  showed  that  the  maximum  strain  exceeded  the  yield 
strain  and  changed  from  the  circumferential  to  the  axial  directions  approximately 
5  to  9  pipe  diameters  ahead  of  the  propagating  crack  tip.  as  shown  in  Rgures  2  and 
3.  These  results  indicated  that  a  long  and  narrow  axial-stretch  zone  preceded  the 
propagating  crack  tip.  This  large  axial  strain,  which  is  much  larger  than  those 
reported  previously,  is  consistent  with  the  observations  of  Shoemaker  et  al.  [9]  and 
Urednicek  [10].  The  narrow  axial-stretch  zone  is  also  the  essential  feature  of  the 
thin  shell  model  of  a  bursting  pipe  by  Freund  et  al.  [11]  who  used  this  model  to 
predict  crack  kinking  at  arrest  due  to  a  mechanical  crack  arrester  [12]. 


Rg.  2.  Strain  rosette  readings 
pressurized  steel  pipe 
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Fig.  3.  Variation  in  principal  strain  direction  derived  from  data  of  Fig.  2 


The  paint,  which  was  sprayed  as  a  background  to  highlight  the  painted  grid  lines, 
acted  as  a  brittle  coating  in  this  test.  The  presence  of  large  axial  strains 
preceding  the  propagating  crack  tip  and  the  alignment  of  the  principal  strain 
direction  with  the  axial  direction  of  the  crack  in  the  region  of  the  crack  tip  was 
graphically  demonstrated  by  the  brittle  coating  crack  pattern  in  Figure  4  where  the 
crack  arrested  in  the  ungrooved  section  of  the  pipe. 

CRACK  TIP 


Figure  4.  Brittle  coating  cracking  pattern  along  crack  path 
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Figure  .4  shows  the  axial  principal  strain  orientation  and  the  region  of  large 
straining  which  extended  on  both  sides  of  the  crack  path  and  at  least  one  inch 
ahead  of  the  crack  tip. 

The  strain  at  the  loci  of  the  crack  ends  of  the  brittle  coating  pattern,  i.e.,  along  the 
isoentatics,  is  estimated  to  be  larger  than  two  percent.  Also  to  be  noted  is  the 
change  in  the  direction  of  the  cracking  pattern  approximately  1/2  in.  away  from 
the  axial  crack  path  in  the  upstream  region.  The  30**  change  in  crack  pattern 
indicates  that  the  principal  strain  direction  abruptly  rotated  upstream  of  the  crack 
tip  with  the  opening  and  stretching  of  the  flaps  and  is  consistent  with  Freund's 
shell  model  of  axial  pipe  fracture  [11]. 

The  width  of  the  brittle  coating  crack  pattern  and  the  lengthening  of  the  painted 
grid  lines  indicated  that  a  large  axial  plastic  straining  occurred  in  a  narrow  region 
surrounding  the  crack  tip.  The  brittle  coating  pattern  thus  provided  a  qualitative 
confirmation  of  the  strain  gage  measurements  and  its  direction  as  well  as  the 
axial  stretch  model  of  the  pipe  flaps. 

The  CTOA  of  a  crack,  which  ran  out  of  the  groove  and  arrested  in  the  ungrooved 
section,  is  almost  triple  that  of  a  running  crack.  This  large  CTOA  suggests  that 
the  crack  flap  motion  continued  unabated  for  a  short  duration  after  the  crack  had 
slowed  to  an  arrest  _ 

tJvnamic  Crack  Curvino  and  Branching 

The  function  of  a  tear  strap  in  a  pressurized  fuselage  is  not  so  much  to  arrest  a 
propagating  axial  crack  by  reducing  the  circumferential  stress  in  the  crack  path, 
but  to  deflect  the  crack  in  the  circumferential  direction.  A  larger  opening 
generated  by  the  deflected  crack  will  then  result  in  a  controlled  depressurization 
of  the  fuselage  and  hence  eliminate  the  driving  force. 

The  mechanics  of  elastic  crack  curving  as  well  as  crack  branching  was  studied  by 
one  of  the  authors  and  his  colleague  under  an  ONR  contract  [13,  14].  The  dynamic 
crack  curving  criterion,  which  is  a  dynamic  extension  of  the  static  criterion  by 
Streit  and  Rnnie  [15],  postulates  that  the  micro-cracks  ahead  of  a  crack  tip 
dictates  the  direction  of  crack  propagation.  The  crack  curving  criterion  assumes 
that  when  the  circumferentiai  stress  within  a  prescribed  crack  tip  region  attains  a 
maximum  value  off  the  axis  of  a  self-similar  crack  extension,  crack  curving  will 
occur.  This  maximum  condition,  which  is  based  on  linear  elastic  fracture 
mechanics  (LEFM),  results  in  a  characteristic  crack  tip  distance,  ro.  in  which  the 
propagating  crack  will  deviate  from  its  axis  by  an  angle  of  9c-  For  a  static  crack, 
these  values  are 
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and  the  crack  curving  angle,  9c,  is 


9c  =  cos-1 


where  Ki  and  Oox  are  the  mode  I  stress  intensity  factor  and  the  remote  stress 
component  or  the  first  nonsingular  term  in  the  elastic  crack  tip  stress  field, 
respectively.  The  elastic  crack  curving  criterion  requires  that  ro  <  rc  for  the  crack 
to  curve  an  angle,  9c,  away  from  its  axis  where  rc  is  a  material  constant  which 
specifies  the  characteristic  crack  tip  region  in  which  the  off  axis  micro-cracks 
are  triggered  and  connected  to  the  main  crack  tip. 


In  the  presence  of  a  large  driving  force,  or  a  large  dynamic  stress  intensity  factor, 
the  crack  will  bifurcate  in  order  to  shed  the  excess  driving  force  and  thus  results 
in  crack  branching  while  the  crack  branching  angle  is  governed  by  the  crack  curving 
criterion.  This  crack  branching  criterion  was  used  successfully  to  correlated  the 
predicted  and  measured  crack  branching  angle,  20c,  and  the  estimated  crack 
branching  stress  intensity  factor,  Kis,  in  a  bursting  steel  pipe  [16) 

Finite  Element  Modeling  of  Axial  Rupture  of  Aircraft  Fuselage 


Axial  rupture  of  an  aircraft  fuselage  differs  significantly  with  the  rupture  of  a 
line  pipe  in  that  the  ratio  of  the  half  crack  length,  a,  to  the  fuselage  radius,  R.  is 
one  to  two  orders  of  magnitude  smaller  than  that  of  the  line  pipe.  Also  the  crack 
flaps  of  the  fuselage  are  initially  loaded  statically  by  the  insulation,  which  acts 
as  a  bladder  to  maintain  the  cabin  pressure  on  the  flaps,  and  later  by  the 
aerodynamic  windage  as  the  flaps  peel  open.  Figure  5  shows  the  axial  and 
circumferential  strain  variations  with  crack  extension  from  the  static  tests 
conducted  in  the  DOT  pipe  rupture  studies  where  the  pipe  pressure  was  maintained 
by  a  neoplene  bladder  reinforced  by  a  Kevlar  insert.  The  small  flap  openings  in 
these  tests,  together  with  static  pressure  acting  on  the  flaps  simulated  the  initial 
phase  of  an  axial  rupture  of  a  fuselage.  As  shown  in  Figure  5,  the  deforming  flap  in 
the  static  line  pipe  test  also  generated  an  axial  strain  which  exceeds  the 
circumferential  strain  at  an  distance  approximately  one  half  of  the  pipe  radius 
ahead  of  the  crack  tip.  This  recurring  theme  in  both  the  static  and  dynamic  axial 
crack  studies  of  the  line  pipe  is  the  large  axial  strain  which  precedes  the  crack 
tip.  This  axial  strain,  which  is  three  to  five  times  larger  than  the  circumferential 
strain,  is  the  cause  of  crack  curving  once  the  propagating  crack  runs  out  of  the 
groove  in  the  line  pipe. 
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In  order  to  estimate  the  axial  strain  ahead  of  the  crack  tip  in  a  rupturing  aircraft, 
a  finite  element  model  of  a  two  bay  crack  in  an  idealized  fuselage,  as  shown  in 
Figure  6,  was  studied.  This  quasi-static  analysis  accounts  for  the  deforming  flap 
geometry  but  ignores  the  stress  wave  effects.  The  quasi-static  analysis  was 
justified  since  the  measure  crack  velocities  in  the  DOT  pipe  rupture  studies  [7,8], 
full-scale  steel  pipe  rupture  studies  [17]  and  dynamic  fracture  testing  of  single¬ 
edge  notched,  7075-T6  and  7178-T6  aluminum  specimens  [18]  all  recorded  crack 
velocities  less  than  ten  percent  of  the  dilatational  wave  velocities. 


Table  1  shows  some  of  the  structural  details  and  the  material  properties  of  the 
idealized  fuselage  considered  in  this  study.  The  fuselage  was  assumed  to  be  a 
pressurized  cylindrical  tube  with  an  uniform  axial  tensile  stress  of  pR/2t.  While 
structural  details,  such  as  butt  and  overlap  joints  were  ignored,  the  frame 
stiffness  was  incorporated  by  the  beam  elements  as  shown  in  Figure  6.  Windage 
loading  on  the  crack  flaps  was  not  considered  due  to  the  lack  of  measured  pressure 
distribution  on  a  fuselage  flap  of  a  moving  airplane. 


Parameter 

Symbol 

Value  1 

Unit 

Material 

2024-T3  Clad  Sheet 

Skin 

Thickness 

t 

0.036 

Inches 

and 

Ultimate  Strength 

Ftu 

62 

ksi 

Yield  Strength 

Fty 

47 

ksi 

Tear 

Young's  Modulus 

E 

1 0.5x1 0^ 

ksi 

Strap 

Fracture  Toughness 

Kc 

91  * 

ksi  Vin 

Fuselage  Diameter 

D 

148 

Inches 

Tear  Strap  Spacing 

S 

1  0 

Inches 

Tear  Strap  Width 

W 

2 

Inches 

Frame  Spacing 

b 

20 

Inches 

Crack  Length 

2a 

40 

Inches 

Limit  Pressure  Differential 

Pllmit 

7.5 

psi 

Fail-Safe  Pressure  Differential 

Pf» 

9.0 

psi 

*  Ref.  MCIC-HB-OIR,  "Damage  Tolerant  Design  Handbook* 


Table  1.  Structural  and  material  property  data  for  aircraft  fuselage  analysis 


Elastic  and  elastic-plastic  static  analyses,  both  based  on  large  deformation  shell 
theory,  were  conducted.  Both  analyses  were  conducted  for  a  stationary  axial  crack 
of  half  crack  length  of  a  a  is  to  20  Inches  at  one  inch  intervals.  As  such,  this 
analysis  does  not  account  for  the  history  of  crack  extension  where  the  crack  would 
grow  from  an  initial  crack  length  of  the  order  of  1  in.  to  a  «  15  -  20  inches.  Thus 
the  effect  of  the  residual  stresses  in  the  unloaded  crack  flaps  is  also  ignored.  In 
all  analyses,  the  center  frame,  as  shown  in  Figure  7  was  assumed  to  be  intact  and 
the  riveted  skin  was  assumed  to  have  torn  off  10  in.  away  from  the  crack.  The 
results  of  this  analysis  is  discussed  in  the  following. 
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TEAR  STRAP 


Rg.  7.  Axial  rupture  of  an  idealized  fuselage 


Elastic  .Analysis 

As  the  bursting  pipe  study  showed,  large  plastic  deformation  precedes  the 
propagating  crack  tip  which  is  followed  by  extensive  plastic  unloading  along  the 
edges  of  the  flaps.  Although  the  elastic  analysis  in  this  section  does  not  model 
the  actual  flap  deformation  process  accurately,  this  zeroth  order  approximation 
provides  an  insight  to  the  problem  at  hand. 

Figure  8  shows  the  elastic  stress  distribution  in  a  cracked  fuselage,  which  is 
pressurized  to  a  gage  pressure  of  7.5  psi  and  an  axial  tensile  stress  of  7.7  ksi, 
with  a  crack  length  of  a  «  16  in.  The  numbers  with  an  asterisk  are  the  computed 
axial  stresses  in  each  finite  element.  The  lower  numbers  represent  the  axial  to 
circumferential  stress  ratios.  The  coarse  finite  element  mesh  of  Figure  6  was 
calibrated  against  a  known  solution,  i.e.  a  central  crack  in  an  infinite  plate  under 
uniaxial  tension  in  this  case,  such  that  the  stress  intensity  factor  and  the  remote 
stress  component,  oqx.  could  be  extracted  from  the  circumferential  and  axial 
stresses  at  the  center  of  the  forward  crack  tip  element.  This  calibration  scheme 
was  essential  for  extracting  meaningful  fracture  data  from  the  coarse  grid  finite 
element  analysis. 
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Rg.  8.  Elastic  stress  distribution  in  a  cracked  pressurized  fuselage 

Table  2  shows  the  mode  I  stress  intensity  factor,  K,  and  the  characteristic 
distance,  ro,  for  half  crack  lengths  ranging  from  a  *  15  through  20  in.  where  the 
crack  was  assumed  to  have  penetrated  halfway  through  the  tear  strap.  Also  shown 
is  the  crack  curving  angle,  dc.  which  was  predicted  by  using  the  crack  curving 
criterion  as  described  in  the  previous  section,  for  these  crack  lengths.  The 
critical  material  property  for  crack  curving,  rc.  was  assumed  to  be  1.0  ~  1.5  mm 
following  similar  analysis  for  steel  pipes  [16].  Note  that  the  high  axial  stress  for 


(deg) 


K 

ihaWi 

144. 

149, 

151. 

157. 

110. 


73.0 

88.3 
109. 
130. 

52.3 


much  of  the  crack  length,  a,  would  have  caused  the  crack  to  curve  if  the  fuselage 
was  not  weakened  by  stress  concentrations  and  multiple  site  damages  along  its 
path.  Thus.  sim'’ur  to  the  2-in.  pipe  rupture  tests,  the  axial  crack  in  the  fuselage 
propagated  straight  ahead  only  in  the  presence  of  the  MSD’s  which  acted  as  a  crack 
guide. 

Exception  to  the  tendency  for  crack  curving  can  be  seen  in  longest  crack  which  was 
assumed  to  have  penetrated  into  the  tear  strap.  In  this  case,  ro  >  rc  and  the 
necessary  condition  for  crack  curving  is  violated,  thus  suggesting  that  the  crack 
will  propagate  axially  once  it  has  fractured  the  tear  strap. 


Figures  9,  10  and  11  show  the  axial  stresses  and  the  axial-to-circumferential 
stress  ratios  in  the  finite  elements  surrounding  the  crack  tip  for  half  crack  length 
of  a  -  15,  18  and  20  in.  Figures  9  and  10  represent  the  states  of  crack  tip  stress 
as  the  crack  approaches  the  tear  strap  and  Figure  11  represents  the  state  of  stress 
when  the  crack  penetrated  the  tear  strap.  Note  that  the  axial*to-circumferential 
stress  ratio  exceeds  unity  for  a  »  18  in.  and  is  a  low  0.6  when  the  crack  had 
penetrated  the  tear  strap. 


The  crack  curving  and  crack  branching  criteria  described  above  obviously  is  not 
applicable  to  the  elastic-plastic  states  of  stress  of  Rgures  9,  10  and  11.  if, 
however,  the  state  of  stress  governs  the  fracture  process,  then  the  crack  curving 
and  crack  branching  criteria  can  be  used  to  estimate  the  onset  of  crack  curving  as 
well  as  the  crack  curving  angle.  The  rational  for  using  such  crude  approximation  is 
analogous  to  the  use  of  plasticity  corrected  stress  intensity  factor  in  LEFM. 


■ 

Fig.  9. 


B  INDICATES  PLASTIC  REGION 
*  INDICATES  AXUL  STRESS  a,  (KSI) 
NUMBERS  INDICATB  ■  -  O,  /  ^ 

WHERE  Oy  CIRCUMFERENTIAL  STRESS 


FRipe 


Elastic-plastic  stress  distribution  in  a  cracked  pressurized  fuselage 
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Fig.  10.  ElastiC'plastic  stress  distribution  in  a  cracked  pressurized  fuselage 


■ 


■  INDICATES  PLASTIC  REQION 
*  INDICATSS  AXIAL  STRESS  O,  (KSI) 
NUMiERS  INDICATE  ■  .  O  jo 

___  *  '' 
WHO*  Oy  CIRCUMFERENTIAL  STRESS 

FRAME 


Fig.  11.  Elastic-plastic  stress  distribution  in  a  cracked  pressurized  fuselage 
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Table  3  shows  the  plasticity  corrected  mode  I  stress  intensity  factor,  K,  the 
remote  stress  component,  oox.  the  characteristic  distance,  ro,  and  the  crack 
curving  angle,  dc  for  two  fuselage  pressures  of  p  =  4.5  and  7.5  psi.  The  results 
show  that  crack  curving  is  inevitable  in  the  skin  before  the  crack  arrived  at  the 
tear  strap. 


p 

(DSil 

a 

(ini 

K 

<T0X 

fksil 

To 

^mm) 

0  c 

re 

=1mm 

(deg) 

rc 

=1 .5mm 

15 

75.5 

62 

0.8 

26.7 

46.5 

16 

78.1 

62 

0.9 

20.9 

44.0 

4.5 

17 

79.2 

62 

0.9 

17.9 

42.9 

18 

73.0 

62 

0.8 

31.2 

48.7 

19 

85.3 

28.4 

5.1 

- 

- 

20 

83.9 

39.1 

2.6 

. 

. 

1  5 

86.2 

61.7 

1.1 

- 

34.9 

16 

88.3  ! 

62 

1.1 

- 

32.7 

7.5 

17 

85.2 

62 

1.1 

- 

36.6 

18 

78.5 

62 

0.9 

19.9 

43.6 

19 

90.2 

24.7 

7.6 

- 

- 

20 

90.8 

36.6 

3.5 

- 

- 

Table  3.  Results  of  elastic-plastic  analysis 


The  fuselage  in  this  analysis  is  a  shell  of  uniform  thickness  without  rivets.  The 
stress  concentrations  at  the  row  of  rivet  holes  in  the  actual  fuselage  will  locally 
increase  the  circumferential  stress  and  simulate  the  groove  in  the  line-pipe 
experiment.  Thus  an  axial  crack  in  the  actuai  fuseiage  wiil  continue  to  propagate 
axially,  particularly  in  the  presence  of  a  multiple  site  damage.  As  the  axial  crack 
in  such  damaged  fuselage  approaches  the  tear  strap,  however,  crack  curving  is 
inevitable  if  the  integrity  of  the  tear  strap  is  maintained  as  shown  by  Figure  10. 
If  the  axial  crack  penetrates  the  tear  strap,  as  shown  in  Figure  11,  then  it  will 
continue  to  propagate  axially  without  curving  through  the  tear  strap. 


Conclusions 

Large  deformation  elastic-plastic  finite  element  analysis  of  a  pressurized 
fuselage  showed  that  a  large  axial  stress  component  exists  ahead  of  an  axial  crack. 
The  computed  crack  tip  stress,  when  combined  with  a  previously  developed  crack 
curving  criterion,  predicted  crack  curving  throughout  much  of  the  skin. 
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Discussions 

The  numerical  analysis  and  results  presented  in  this  paper  should  be  considered 
a  first  order  approximation  to  the  complex  state  of  stress  in  an  actual 
fuselage.  The  analysis  could  be  easily  upgraded  by  incorporating  the  history 
effect  of  plastic  deformation  by  incrementally  advancing  the  axial  crack  tip. 

A  more  difficult  task  is  to  develop  a  crack  curving  and  crack  branching  criteria 
in  the  presence  of  large  plastic  deformation.  If,  for  example,  a  maximum 
strain  criterion  is  used  for  crack  curving,  then  crack  curving  will  not  occur 
until  the  crack  tip  is  close  to  the  tear  strap  as  indicated  by  the  axial-to- 
circumferential  stress  ratios  in  Figures  9.  10  and  11. 
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